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Abstract 

Measurements of cross sections for beauty and charm events with dijets and a muon in 
the photoproduction regime at HERA are presented. The data were collected with the 
HI detector and correspond to an integrated luminosity of 179 pb _1 . Events with dijets 
of transverse momentum P^ 11 > 7 GeV and P^ a > 6 GeV in the pseudorapidity range 
— 1.5 < rf et < 2.5 in the laboratory frame are selected in the kinematic region of photon 
virtuality Q 2 < 2.5 GeV 2 and inelasticity 0.2 < y < 0.8. One of the two selected jets must 
be associated to a muon with Pj,> 2.5 GeV in the pseudorapidity range —1.3 < rf< 1.5. 
The fractions of beauty and charm events are determined using the impact parameters of the 
muon tracks with respect to the primary vertex and their transverse momentum relative to 
the axis of the associated jet. Both variables are reconstructed using the HI vertex detector. 
The measurements are in agreement with QCD predictions at leading and next- to-leading 
order. 
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1 Introduction 



The photoproduction of heavy quarks at the HERA electron-proton collider is of particular inter- 
est for testing calculations in the framework of perturbative quantum chromodynamics (pQCD). 
The masses m b and m c of the beauty and charm quarks, as well as the transverse momentum 
of a jet, provide a hard scale, which is necessary for the calculation of pQCD predictions 0]. 
The dominant process for beauty and charm production at HERA is boson-gluon fusion (BGF): 
7<? QQX, with Q — b,c. In the kinematic regime of photoproduction, where the photon 
virtuality Q 2 is small, two classes of processes contribute to BGF. In direct photon processes, 
the photon emitted from the electron^ enters the hard BGF process directly. In resolved pho- 
ton processes, the photon fluctuates into a hadronic state before the hard interaction and acts 
as a source of partons, one of which takes part in the hard interaction. At leading order (LO) 
pQCD resolved photon processes are expected to contribute significantly in the photoproduction 
region. 

Heavy quark photoproduction has been investigated at HERA using different analysis tech- 
niques in different regions of phase space. The lifetime or mass of heavy flavoured hadrons 
[]2]-(71, semileptonic decays JMMHl or the full reconstruction of a D or D* meson ll9l[T0ll are ex- 
ploited to perform the measurements. In general, the measured cross sections agree reasonably 
well with the theoretical predictions. In the HI measurement performed with HERA I data 0, 
where beauty photoproduction was investigated using two jets and a muon in the final state, the 
data were found to be described reasonably well by the pQCD calculations at next-to-leading 
order (NLO), except in the region of low transverse momentum of the muon 2.5 < < 3.3 
GeV and of low leading jet transverse momentum 7.0 < P^ tl < 10.0 GeV, where the NLO 
prediction was lower than the data. Comparable measurements were performed by the ZEUS 
collaboration [|6]I3, covering a similar phase space. Here, a good description by the NLO QCD 
predictions is found, also at low transverse momenta of the leading jet and the muon. 

This paper reports on beauty and charm measurements of cross sections for photoproduc- 
tion of events with two jets and a muon, where the muon is associated with one of the jets. The 
analysis exploits the lifetime and the mass of heavy flavoured hadrons as in the former HERA 

1 beauty production analysis @. The measurement is additionally performed for charm photo- 
production. Increased statistics due to increased luminosity and extended phase space, a better 
understanding and description of the HI vertex detector [fill , which is crucial for this analysis, 
as well as reduced systematic uncertainties compared to the previous HI analysis make a more 
detailed test of pQCD predictions possible. Going beyond the earlier HERA I analysis, the 
cross sections as a function of the azimuthal angular difference A</4> between the two leading 
jets, which are sensitive to higher order corrections, are also measured. Furthermore, cross 
sections are measured in two different regions of phase space, which are either enriched by 
resolved or direct photon processes. 

2 QCD Calculations 

The data presented here are compared with LO calculations supplemented by parton show- 
ers as well as with NLO calculations. The calculations are performed using either collinear 

'in this paper the term 'electron' is used generically to refer to both electrons and positrons. 
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factorisation, which is based on the DGLAP evolution equations lfT2ll . or the k t - factorisation 
approach, which employs the CCFM [13] evolution equations. In the collinear approach trans- 
verse momenta obtained through initial state QCD evolution are neglected and all the transverse 
momenta are generated in the hard scattering process, i.e. the partons entering the hard inter- 
action are collinear with the proton. Effects from non-zero transverse momenta of the gluons 
enter only at NLO. In the /^-factorisation ansatz the transverse momenta of gluons k t entering 
the hard interaction are already included at leading order both in the off- shell matrix element 
and the k t - dependent unintegrated gluon density [fT4l . Therefore, corrections appearing only 
at higher order in collinear factorisation are partially included at LO in the A; r factorisation ap- 
proach. 

For beauty and charm photoproduction two classes of processes occur, the direct photon 
processes and the resolved photon processes. The distinction between these two classes depends 
on the factorisation scheme and the order in which the calculation is performed. 

The production of heavy quarks is calculated either in the massive scheme, where heavy 
quarks are produced only perturbatively via boson gluon fusion, or in the massless scheme, 
where heavy quarks are treated as massless partons. These two schemes are expected to be 
appropriate in different regions of phase space [15]: the massive scheme is expected to be 
reliable when the transverse momentum P T of the heavy quarks is of similar size compared to 
the heavy quark mass mg, whereas the massless scheme is expected to be valid for P T ^> mq. 

An overview of the parameters used in the Monte Carlo (MC) generators is given in tabled] 
The following MC generators are used: 

Pythia: The MC program Pythia 6.4 O is based on LO QCD matrix elements with 
leading-log parton showers in the collinear factorisation approach. Pythia includes 
both direct photon gluon fusion and resolved photon processes. In the resolved pho- 
ton processes either a beauty or a charm quark or a gluon from the photon enters the 
hard scattering. In the inclusive mode of Pythia used here beauty and charm quarks are 
treated as massless partons in all steps of the calculation in both types of processes. Three 
samples are generated containing photoproduction events for the processes ep — > bfiX, 
ep — >■ cjj,X and ep — > qjX where q is a light quark of flavour u, d or s and j denotes 
a jet. The latter sample is generated without specifically requiring a muon in order to 
use it for studying the background arising from muon candidates which originate from 
sources other than beauty or charm processes. The hadronisation process is simulated 
according to the Lund string model [flD . For the systematic uncertainty arising from the 
fragmentation model, additional samples are generated using the Peterson fragmentation 
function [18] for heavy quarks. 

Cascade: The MC program Cascade 2.0 [[191 is used for simulating events based on LO 
QCD calculations in the k t -factorisation approach. Only the direct boson gluon fusion 
process is implemented using off-shell matrix elements. Higher order QCD corrections 
are simulated with initial state parton showers applying the CCFM evolution [|T3l . Here, 
two samples containing the processes ep — > bfxX and ep — > c[iX are generated. The 
unintegrated PDFs of the proton from set AO [20] are used. The hadronisation of partons 
is performed with the Lund string model as implemented in Pythia. 
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Pythia 


Cascade 


Herwig 


MC@NLO 


Version 
Evolution scheme 


6.4 
DGLAP 


2.0 
CCFM 


6.510 
DGLAP 


Herwig 6.510 
DGLAP 


m b [GeV] 
m c [GeV] 


4.75 
1.5 


4.75 
1.5 


4.75 
1.5 


4.75 
1.5 


Proton PDF 
Photon PDF 


CTEQ6L1 [26] 
SaS2D [27] 


AO 1201 


HERAPDF 1.0 [28] 
GRV-G LO [29] 


HERAPDF 1.0 
GRV-G LO 


Renorm. scale pL R 
Factor, scale /xp 




V § + Qi 

Hr 


V>R 




Fragmentation 


Lund 


Lund 


cluster 


cluster 




a = 0.437 
b = 0.850 
Peterson 
e fe = 0.006 
e c = 0.06 


a = 0.437 
b = 0.850 







Table 1: Parameters used in the QCD calculations Pythia, Cascade, Herwig, and 
MC@NLO. Here, Pt stands for the transverse momentum, mi, rri2 are the masses of the two 
outgoing partons from the hard process, and tuq denotes the heavy quark mass. The variable s 
represents the invariant mass of the QQ system and Q± stands for its transverse momentum. 



Herwig: The MC program Herwig 6.510 [21] is used to simulate events based on collinear 
factorisation and massless LO QCD calculations. Herwig includes both direct photon 
gluon fusion and resolved photon processes. The hadronisation of partons is performed 
using the cluster fragmentation [1221 . 

MC@NLO: The MC@NLO program for HERA [1] combines a collinear NLO calculation of 
heavy flavour production in the massive approach [|23l with parton showers and hadro- 
nisation as described in II241 . The direct and resolved parts of the cross section are cal- 
culated separately. MC@NLO uses parton showers applying the DGLAP evolution to 
simulate higher order contributions and the cluster fragmentation [|22l as implemented in 
Herwig (2TJ. The quark masses are set to m c = 1.5 GeV and m b = 4.75 GeV. The 

central value of the renormalisation scale fiR is set to J nig + P T , while the factorisa- 
tion scale /xp is 2/xr. As an estimate of the theoretical uncertainties on the NLO QCD 
predictions the scales jiR and /xp are varied separately by factors of 0.5 and 2, while m b 
is changed by ±0.25 GeV and m c is changed by ±0.2 GeV. The resulting variations are 
added in quadrature to obtain the resulting theoretical uncertainty. 

Pythia and Cascade are used to simulate detector effects in order to determine the ac- 
ceptance and the efficiency and to estimate the systematic uncertainties associated with the 
measurement. The generated events are passed through a detailed simulation of the detector 
response based on the GEANT simulation program [25] and are processed using the same re- 
construction and analysis chain as is used for the data. 
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3 HI Detector 



Only a short description of the HI detector is given here including the most relevant detector 
components for this analysis. A more complete description may be found elsewhere [|30ll3Ti . 
A right-handed coordinate system is employed at HI, with its origin at the nominal interaction 
vertex, its z-axis pointing in the proton beam direction and its x(y) axis pointing in the hori- 
zontal (vertical) direction. Polar (9) and azimuthal ((f)) angles are measured with respect to this 
reference system. The pseudorapidity 77 is related to the polar angle 9 by 77 = — mtan(0/2). 

Charged particles are measured in the central tracking detector (CTD) with a transverse 
momentum resolution of a(P T )/ P T m 0.5%Pt/GcV © 1.5% [|32l . This device consists of 
two cylindrical drift chambers (CJC) interspersed with a drift chamber designed to improve the 
^-coordinate reconstruction. A multiwire proportional chamber mainly used for triggering is 
located in front of the inner CJC. The CTD is operated in a uniform solenoidal 1.16 T magnetic 
field, enabling the momentum measurement of charged particles over the polar angular range 
20° < 9 < 160°. The efficiency for finding tracks in the CTD is greater than 99%. 

The CTD tracks are linked to hits in the vertex detector, the central silicon tracker (CST) 
ll33l . to provide precise spatial track reconstruction. The CST consists of two layers of double- 
sided silicon strip detectors surrounding the beam pipe, covering an angular range of 30° < 9 < 
150° for tracks passing through both layers. The information of the z-coordinate of the CST hits 
is not used in the analysis presented in this paper. For CTD tracks with CST hits in both layers 
the transverse distance of closest approach (DCA) to the nominal vertex in x — y, averaged over 
the azimuthal angle, is measured to have a resolution of 43 fim © 51 /mi/ (P^/GeV), where the 
first term represents the intrinsic resolution (including alignment uncertainty) and the second 
term is the contribution from multiple scattering in the beam pipe and the CST. The efficiency 
for linking hits in both layers of the CST to a CTD track is around 84%. 

The track detectors are surrounded in the forward and central directions (4° < 9 < 154°) by 
a finely grained liquid argon calorimeter (LAr) and in the backward region (153° < 9 < 178°) 
by a lead-scintillating fibre calorimeter (SpaCal) both with electromagnetic and hadronic sec- 
tions. These calorimeters provide energy and angular reconstruction for final state particles 
from the hadronic system. In the LAr electromagnetic shower energies are measured with 
a precision of a(E)/E = 11%/v/P/GeV © 1% and hadronic energies with a(E)/E = 
50%/ 'y/E/GeV © 1%, as determined in test beam measurements. The energy resolution for 
electromagnetic showers in the SpaCal is a(E)/E = 7.1%/ \jEj GeV © 1%, as determined in 
test beam measurements [|34l . 

The calorimeters are surrounded by the muon system. The central muon detector (CMD) 
is integrated in the iron return yoke of the superconducting coil and consists of 64 modules, 
which are grouped in the forward endcap (5° < 9 < 35°), the forward and backward barrel 
(35° < 9 < 130°) and the backward endcap (130° < 9 < 175°). Muon candidates are identified 
by requiring a geometrical matching of a CMD track segment with a CTD track. 

The luminosity determination is based on the measurement of the Bethe-Heitler process 
ep — > ep7, where the photon is detected in a calorimeter located downstream of the interaction 
point in the electron beam direction at z = — 103 m. 
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4 Experimental Method 



The data were collected with the HI detector at the HERA collider during the years 2006 and 
2007 and correspond to an integrated luminosity of C = 179 pb -1 . The beam energies were 
E e = 27.6 GeV and E p = 920 GeV for electrons and protons, respectively, resulting in a centre- 
of-mass energy of ^fs m 320 GeV. The trigger requires a track segment in the muon system and 
track activity in the central jet chamber. A detailed account of the present analysis can be found 
in 11351 . A summary of the kinematic range and the definition of the measurement is given in 
tabled 

4.1 Photoproduction Event Selection 

Events in the photoproduction regime are selected by requiring that no isolated high energy 
electromagnetic cluster, consistent with a signal from a scattered electron, is detected in the 
LAr and SpaCal calorimeters. This limits the photon virtuality to values of Q 2 < 2.5 GeV 2 . 
The inelasticity y is reconstructed using the relation y = ^2 h (E — P z )/2E e ||36ll . Here, the 
sum includes all particles of the hadronic final state (HFS), while E denotes their energies and 
P z stands for the ^-components of their momenta. The HFS particles are reconstructed using a 
combination of tracks and calorimeter deposits in an energy flow algorithm that avoids double 
counting [37J. The inelasticity in this analysis is restricted to 0.2 < y < 0.8. 

4.2 Muon Reconstruction and Selection 

Muon candidates are identified as track segments in the barrel and endcap parts of the instru- 
mented iron. The iron track segments must be well matched to a track reconstructed in the 
CTD. At least two CST hits in the r — <p plane have to be associated with the muon track. The 
combined CTD-CST track in r — is required to have a fit probability of at least 10%. The 
muon momentum is reconstructed using the CTD-CST track information. The CST hit require- 
ments for the muon track restrict the allowed range of ep interaction vertices along the z-axis 
to \z vtx \ ^ 20 cm. Events are selected with at least one muon candidate reconstructed in the 
instrumented iron having a pseudorapidity within —1.3 < ry M < 1.5 and a transverse momen- 
tum of > 2.5 GeV. If more than one muon candidate is found, the one with the highest 
transverse momentum is selected and other candidates are ignored. In 5.4% of the events after 
the full selection more than one muon is found. 



4.3 Jet Reconstruction and Selection 

Jets are reconstructed using the inclusive longitudinally invariant k T algorithm in the massless 
Pt recombination scheme and with the distance parameter Rq — 1 in the 77 — </> plane [38]. The 
algorithm is applied in the laboratory frame using all reconstructed HFS particles including the 
muon candidate. A jet is defined as a /i-jet if the selected muon candidate lies within a cone 
of radius 1 about the jet axis in the rj — (j) plane. The efficiency for this matching amounts to 
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Photoproduction of 6(c) — > fxjjX 


Kinematic range 


Q 2 < 2.5 GeV 2 
0.2 < y < 0.8 


Event selection 


P£> 2.5 GeV 
-1.3 < 77" < 1.5 
pj, etl > 7 GeV 
P^ 2 > 6 GeV 
-1.5 < rf et < 2.5 


Event sample 


N events = 6807 
C = 179 pb" 1 



Table 2: Definition of the kinematic range of the measurement and event yield for the data 
sample collected in the years 2006 and 2007. The variables are measured in the laboratory 
frame. 

about 90% and is consistent for data and MCs. The jet with the highest P T is referred to as 
jetl, while the second highest is called jet2. Events with at least two jets are selected, where 
the leading two jets are required to be in the angular range —1.5 < rf et < 2.5 and to have a 
transverse momentum of P^ tl > 7 GeV and P^ 12 > 6 GeV. One of the two selected jets must 
be classified as a /i-jet. The Monte Carlo simulation is used to define hadron level jets, which 
consist of stable particles including neutrinos, but excluding the scattered electron, before they 
are passed through the simulation of the detector response. 



4.4 Separation of Direct and Resolved Processes 

The fraction of the photon energy entering the hard interaction is estimated using the observable 

^.obs. 

r0bs = E jet iiE-P z ) + E jea (E-P z ) 
Ehfs(E-P z ) 

where the sums in the numerator run over the particles associated with the two jets and the 
one in the denominator over all detected hadronic final state particles. For direct processes x° bs 
approaches unity, as the hadronic final state consists only of the two hard jets selected in the 
present analysis and the proton remnant in the forward region only has a minor contribution to 
J2hfs(E — P z )- In resolved processes x° bs can have smaller values. 



4.5 Flavour Separation 

The flavour of an event is defined as the hadron flavour of the /i-jet. The measured cross sections 
are proportional to the rate of events with a muon and a dijet system rather than the rate of muons 
or jets. 

The separation of b, c and light quark (uds) events is only briefly described here. The 
procedure closely follows that described in [Q. The separation is performed using the properties 
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of the muon track associated to the /U-jet. The impact parameter 5 of a muon track is the 
transverse distance of closest approach (DCA) of the muon track to the beam spot point, which 
is the position of the beam interaction region in x and y. The beam spot is derived from tracks 
with CST hits averaging over many events and is updated regularly to account for drifts during 
beam storage. Muon tracks with 5 > 0.1 cm are rejected to suppress contributions from the 
decays of long-lived strange particles. If the angle a between the azimuthal angle of the yU-jet 
and the line joining the primary vertex to the point of muon closest approach is less than 90°, 
the impact parameter is defined as positive. It is defined as negative otherwise. The transverse 
momentum P^ el of the muon relative to the /i-jet axis is also sensitive to the quark content of 
the event sample and used together with the impact parameter for the flavour separation. 

The fractions of events with beauty, charm and light quarks, ft, f c and fi, are obtained by a 
binned likelihood fit ||39ll in the 5 — P^ el plane. Following 11391 , a likelihood ratio is calculated 
based on Poisson statistics. The fit is performed separately for each individual measurement 
bin, while the total cross sections are determined using the fractions obtained from a fit to the 
complete event sample. The uds (light), c and b Pythia Monte Carlo simulation samples 
are used as templates. Only the statistical errors of the data and the Monte Carlo simulations 
are considered in the fit. As a cross check, all fits are also performed using one-dimensional 
distributions of P^ el and 5 separately. These two one-dimensional fits give a compatible b 
fraction in all measurement bins within the statistical error. The one-dimensional P^ el fit does 
not allow a determination of the charm fraction. 

The two distributions that are used in the flavour separation are shown in figure [Q The 
distribution of the impact parameter 5 shown in figured] (a) is symmetric for uds events while 
b and c events contribute more at large positive values of 5. Therefore, the fit of this variable 
can distinguish the three different quark contributions. As shown in figure \T\ (a), the sum of 
all three fitted contributions in the Monte Carlo simulation is able to describe the data quite 
well. This description is achieved by a better understanding of the detector and an improved 
detector simulation with regard to signal heights, noise levels and dead strips in the CST [fTTTl . 
the inclusion of effects from alignment imperfections, and the description of the dead material 
in front of the CST and CJC. Therefore, a further smearing of measured track parameters in the 
simulation is not necessary, as it was done in the former HI analysis using HERA I data 0. 
The P^ el distribution is shown in figure [11(b). The uds and c distributions are very similar and 
peak at low values of P^ el , while the b events contribute more at higher values of P^ el . The sum 
of all quark contributions in the Monte Carlo simulation is able to describe the data reasonably 
well. 

The fitted parameters /&, f c and /; for the whole kinematic range are: 

f b = (26.0 ± 1.2)%, 
f c = (48.6 ± 2.5)%, 
ft = (25.3 ± 2.6)%. 

The x 2 /ndf is found to be 0.76 for the total sample. 

Control distributions for the data sample in comparison to the Monte Carlo simulations 
are shown in figure |2] All selection cuts are applied. The data are compared with the MC 
contributions from beauty, charm and light quark events and their sum with the relative fractions 
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taken from the fit as discussed above. The number of events in the simulation is normalised to 
the one of the data. It is observed that the shapes of the MC contributions are rather similar for 
beauty, charm and light quark events for the distributions shown here. For the determination 
of the detector corrections a reweighting of Acfiti and the transverse momentum of the leading 
selected track in the event is performed on hadron level in the Monte Carlo simulation to provide 
a better description of the data. Only small deviations between data and MC are observed, such 
as in the forward 77^ region in figure |2](b) and in the very small A<^ region in|2](d). 



4.6 Cross Section Determination 



Total and differential visible beauty and charm ep cross sections are measured in the photopro- 
duction regime. The fitted fractions ft, and f c are converted to cross sections au c ) in each bin 
using 



f at \jMCgen 
_ Jb(c)^Data^ b{c) 

a Kc) ~ rMMCrec " ^ 

Here, N^ata and N^ rec represent the number of data or Monte Carlo simulation events passing 

all selection cuts on reconstruction level. The variable N^' gen stands for the events selected 
in the Monte Carlo simulation on hadron level and C denotes the luminosity of the data. The 
differential cross sections are obtained by dividing by the bin width. 



5 Fake Muon Rate 



All backgrounds to semi-muonic b and c decays are called fake muons here. These contributions 
are modelled using MC simulation and originate mainly from uds events, with a small fraction 
from b and c events. Three sources of fake muons are considered: 



• Hadrons which reach the muon detector and are misidentifled as muons. According to 
the fully inclusive Pythia MC, 1.6% of the selected b events and 2.3% of the selected c 
events originate from hadron misidentification. 

• Muons which do not originate from b or c hadron decays, but from other hadrons such as 
kaons and pions. This background source is denoted as inflight decay in the following. 
According to the full inclusive Pythia MC, 0.9% of the b and 0.7% of the selected c 
events originate from inflight decays. 

• Cosmic ray muons which coincide in time with real ep events. About 1% of the selected 
muon candidates are rejected as cosmic ray muons based on timing information from the 
CTD [35]. The remaining background from cosmic ray muons is negligible. 
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The probability to fake a muon depends on the particle species. The kinematic distributions of 
different particle species differ, but the uds MC is used as one single template. Therefore the 
fake probabilities for the most important particle species are studied in data and MC. The fake 
muon probabilities from misidentification and inflight decays P£ for hadron h are defined as: 

,. # fake muons 
h # all hadrons ' 

The fake muon probability is investigated for pions originating from Kg — > n + n~ decays and 
amounts to Pjf ~ 0.05 in the data. For protons the decay channels A — > pn~ and A — >■ p7r + 
are used and the fake muon probability is found to be, fa 0.04 in the data. For K ± mesons 
from the decay D* ± — > D°7rf low — >■ (K T n ± )7if low a fake muon probability of P£ ± « 0.01 
is measured in the data. It is observed that the pion and proton fake muon probabilities in 
the data are not described by the Monte Carlo simulation. They are reweighted in the Monte 
Carlo simulation by factors of 2.0 and 1.9, respectively, to match the data. The K ± fake muon 
probability in the data and the simulation are in agreement. The misidentified muon events 
remain in the event sample. 



6 Systematic Uncertainties 

The following uncertainties are taken into account in order to evaluate the systematic errors. 

• The trigger efficiencies are determined using independent trigger channels in the DIS 
regime since no independent triggers exist in photoproduction. The uncertainty is esti- 
mated by the difference between the efficiency found in the data and the simulation. The 
Monte Carlo simulation is reweighted to match the data trigger efficiency. 

• The efficiency for the identification of the muons is determined using a high statistics 
sample of events of elastically produced J/tp mesons ||40~1 . The efficiency is known to a 
precision of 4%. 

• The track efficiency of the CTD is known to ±1% and that of the CST to ±2%. The 
uncertainty due to the track efficiencies is estimated by varying the efficiencies of the 
CTD and CST correspondingly. 

• The integrated luminosity is known to a precision of 4%. 

• The uncertainty due to the resolution of the impact parameter 5 of the muon tracks is 
estimated by varying the resolution by an amount that encompasses any difference be- 
tween the data and the simulation. This is achieved by applying an additional Gaussian 
smearing in the Monte Carlo simulation of 200 /im to 5% of randomly selected tracks and 
12 fim to the rest. 

• The uncertainty on the cross section arising from the uncertainty on the reconstruction of 
4>j et is estimated by shifting its value by ±2°. 
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Systematic error source 


Beauty Aer/er[%] 


Charm Aa/a[%} 


Trigger efficiency 


4 


4 


Muon identification 


4 


4 


Track finding efficiency 


3 


3 


Luminosity 


4 


4 


5 Resolution 


3 


2 


Jet axis 


4 


2 


Hadronic energy scale 


3 


5 


Physics model 


3 


1 


Fragmentation 


3 


4 


Fake muon background 


1 


1 


Total 


10.5 


10.4 



Table 3: Summary of the systematic uncertainties of the beauty and charm cross sections. 

• The uncertainty arising from the hadronic energy scale is estimated by changing it by 
±1% for the complete hadronic final state. 

• The dependence of the measurement on the physics model used for the templates repre- 
senting different QCD evolution schemes is estimated by replacing the Pythia b and c 
Monte Carlo templates with CASCADE. 

• The uncertainty on the cross section arising from the uncertainty of the parton fragmen- 
tation model is estimated by replacing the Pythia Monte Carlo samples using the Lund 
fragmentation function with samples based on the Peterson fragmentation function. 

• The uncertainty arising from fake muon background is estimated by not applying the 
weights that have been found in the fake rate probabilities for K and A decays. 

• The impact of the reweighting on the cross sections are investigated and found to be 
negligible. 

The effect of the listed experimental uncertainties are estimated by varying the relevant 
variables in the Monte Carlo simulation or by modifying the corresponding efficiencies in the 
cross section calculation. The difference between the obtained cross sections with and without 
the change result in the measurement systematic uncertainties, which are summarised in table 
[3l The individual effects of the above experimental uncertainties are combined in quadrature, 
yielding a total uncertainty of 10.5% and 10.4% on the measured b and c cross sections, respec- 
tively. The systematic uncertainties as derived from the integrated sample are applied to each 
analysis bin in order to avoid statistical fluctuations. The largest contribution to this uncertainty 
for the c measurement arises from systematics attributed to the hadronic energy scale (5%). The 
systematic errors of the b analysis are not dominated by a single source. 
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0"vis(ep — > ebbX — » ejj/iX') [pb] 


cr vis (ep -> eccX ->■ ejj/iX') [pb] 


HI Data 


43.3 ±2.1 (stat.) ±4.5 (sys.) 


81.3 ±4.3 (stat.) ±8.5 (sys.) 


Pythia 
Cascade 

Herwig 
MC@NLO 


35.3 
29.0 
20.6 
33-41^ 


94.3 
76.8 
58.5 
58.6i??;i 



Table 4: Total visible measured beauty and charm cross sections along with their statistical and 
systematic errors. The total predictions from Pythia, Cascade, Herwig, and MC@NLO 
are also shown. The MC@NLO predictions are given with their theoretical uncertainties. 

7 Results 

The cross sections for b and c in photoproduction using semi-muonic decays in dijet events are 
measured. The cross sections are determined for the phase space defined by the kinematic range 
and the event selection cuts presented in table |2] The measured cross sections are compared to 
the expectations of the MC programs Pythia, Cascade, Herwig, and MC@NLO. The total 
measured and predicted cross sections are listed in table @] Pythia shows the highest normali- 
sation of the three LO MCs, while the normalisation of Cascade is below the one of Pythia 
and Herwig has the lowest normalisation. For the beauty measurement, the Pythia predic- 
tion is closest to the data and gives the best description of the three LO MC predictions. The 
beauty and charm data cross sections tend to be underestimated by the MC@NLO predictions 
but are in agreement within the errors. The precision of the measured cross sections are much 
higher than the ones of the theory predictions shown here. 

The beauty and charm cross sections are measured differentially as a function of the trans- 
verse momentum of the leading jet P^ n and of the muon P£ , the pseudorapidity of the muon 
?7 M , the momentum fraction x° bs carried by the photon entering the hard interaction and the az- 
imuthal angular difference A<^ between the two leading jets. The measurements are performed 
for the full sample, as well as for direct and resolved enriched processes separately. The distinc- 
tion is performed by the variable x° bs , which leads to enriched resolved processes in the region 
x° bs < 0.75 and direct photon enriched processes for x° bs > 0.75. The cross sections for the 
beauty measurements are given in tables |5]|7] and shown in figures [3][5] 

In the case of beauty production the models provide a good description of the measured 
cross sections in terms of shape in all distributions. For the LO MCs this is true for the full 
sample, as well as for the direct and resolved enriched regions. The cross sections for direct 
enriched processes are well described in shape, but tend to be underestimated by MC@NLO, 
while for resolved enriched processes a reasonable agreement is observed both in shape and 
normalisation. 

In the analysis of semi-muonic b decays in dijet events with HERA I data [5] an excess of 
data compared to the NLO predictions of the FMNR program was observed in the first bin of 
P^ and P^ n . In this analysis the NLO predictions are provided by MC@NLO which is based 
on the FMNR parton level calculations. Also in this analysis, the NLO predictions lie below 
the data in the first bin of Pj. and P^ etl , but they are consistent with the data within 2cr of the 
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experimental and theoretical uncertainty. Whereas in enriched direct processes, the data tend to 
be underestimated by the MC@NLO prediction in the first bin of P£ and the first bin of , 
in enriched resolved processes, no such effect is visible. 

The cross sections as a function A<^ show a significant contribution away from the back- 
to-back configuration at A(^j~ 180°. Such a configuration can be described by models which 
include significant contributions from higher order QCD radiation or a transverse momentum 
of the gluon in the initial state. This distribution is reasonably well described by all models. In 
direct and resolved enriched processes, this observation also holds. 

The measured charm cross sections are presented in tables ISlfTOl and figures [6]|8] The dis- 
tributions are reasonably well described by all models. Similar to recent observations in HI 
measurements of the photoproduction of D* mesons [41 J, the central value of the MC@NLO 
calculations tend to be lower than the measured charm cross sections. 



8 Conclusions 

Beauty and charm photoproduction cross sections for events with dijets and a muon are mea- 
sured using the data collected by the HI detector at HERA. Compared to the previous HI beauty 
measurement [5], the analysis profits from a three times larger luminosity of the data sample, an 
extended phase space as well as improved understanding of the HI vertex detector. The flavour 
composition of the event sample is determined by the transverse momentum of the muon rel- 
ative to the jet axis of its associated jet P^ el and by its impact parameter 5. Total visible and 
differential cross sections are measured, and the results are compared to leading order QCD 
models provided by Pythia, CASCADE, and Herwig, as well as to the next-to-leading order 
calculations provided by MC@NLO. At low values of P£ and P^ 1 , the present beauty mea- 
surement does not show a significant excess as observed by the previous HI measurement with 
respect to the NLO calculation. In general the predictions are in reasonable agreement with the 
beauty and charm measurements. 
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HI Beauty Dijet Muon Cross Sections 

&vis(ep —> ebbX — > ejj/iX) 


P£ range 
[GeV] 


da/dP? stat. sys. 
[pb/GeV] 


/&±stat. 


2.5 3.3 
3.3 4.7 
4.7 15.0 


25.3 2.2 2.8 
12.0 0.9 1.3 
0.772 0.078 0.085 


0.21 ±0.02 
0.28 ± 0.02 
0.31 ±0.03 


rf range 


dajdrf stat. sys. 
[pb] 


/&±stat. 


-1.3 -0.3 
-0.3 0.0 
0.0 0.3 
U.o U.O 
0.6 1.5 


9.9 1.2 1.1 

19.6 2.2 2.2 

21.7 2.3 2.4 

OQ A OA Of? 

Z6A ZA Z.b 
14.3 1.4 1.6 


0.23 ±0.03 
0.24 ± 0.03 
0.25 ±0.03 

n on _i_ n no 

0.26 ±0.03 


P^ etl range 
[GeV] 


da/dP^ tl stat. sys. 
[pb/GeV] 


/ 6 ±stat. 


7 11 
11 lo 

15 38 


4.53 0.41 0.50 
o.4U u.Zl u.ol 
0.469 0.038 0.052 


0.24 ±0.02 

n occ _i_ n no 
U.Zo ± U.Uz 

0.30 ±0.02 


A<^> range 
[deg] 


da/dA(f)ti stat. sys. 
[pb/deg] 


/&±stat. 


155 
155 173 
173 180 


0.0576 0.0063 0.0061 
1.01 0.07 0.11 
2.17 0.17 0.24 


0.23 ±0.03 
0.26 ±0.02 
0.26 ±0.02 


x° bs range 


da/dx ^ stat. sys. 
[pb] 


/ 6 ±stat. 


0.0 0.4 
0.4 0.75 
0.75 1.0 


10.6 2.7 1.2 
35.6 3.4 3.9 
103.5 5.9 11.4 


0.17 ±0.04 
0.23 ±0.02 
0.29 ±0.02 



Table 5: Bin averaged differential cross sections for beauty photoproduction of dijet events 
using semi-muonic decays in bins of P^, rf, P^ , A</4>, and x° bs with their statistical and 
systematic uncertainties. The fit parameter f b is given including its statistical error. 
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HI Beauty Dijet Muon Cross Sections, a;° bs > 0.75 

a vis(^P ebbX — > ejjjiX) 


P% range 
[GeV] 


da/dP? stat. sys. 
[pb/GeV] 


/&±stat. 


2.5 3.3 
3.3 4.7 
4.7 15.0 


14.8 1.5 1.6 
6.7 0.6 0.7 
0.487 0.068 0.054 


0.24 ±0.02 
0.28 ±0.03 
0.31 ±0.04 


rf range 


dajdrf stat. sys. 
[pb] 


/&±stat. 


-1.3 -0.3 
-0.3 0.0 
0.0 0.3 
0.3 0.6 
0.6 1.5 


6.5 0.9 0.7 

11.3 1.6 1.2 

14.4 1.8 1.6 
13.9 1.7 1.5 

6.1 1.0 0.7 


0.24 ±0.03 
0.24 ±0.03 
0.29 ±0.03 
0.32 ±0.04 
0.25 ±0.04 


P|? range 
[GeV] 


dajdP^ stat. sys. 
[pb/GeV] 


/ 6 ±stat. 


7 11 
11 15 

15 38 


1.81 0.24 0.20 
2.05 0.21 0.23 
0.361 0.032 0.040 


0.21 ±0.03 
0.26 ±0.03 
0.35 ±0.03 


Acft range 
[deg] 


da/dAcft stat. sys. 
[pb/deg] 


/ 6 ±stat. 


155 
155 173 
173 180 


0.0084 0.0030 0.0016 
0.617 0.052 0.068 
1.640 0.144 0.180 


0.18 ±0.04 
0.28 ±0.02 
0.30 ±0.03 



Table 6: Bin averaged differential cross sections for beauty photoproduction of dijet events 
using semi-muonic decays for x° bs > 0.75 in bins of P£, 77^, P^ etl , and A^j with their statistical 
and systematic uncertainties. The fit parameter f b is given including its statistical error. 
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HI Beauty Dijet Muon Cross Sections, a;° bs < 0.75 

a vis(^P ebbX — > ejjjiX) 


P? range 
[GeV] 


da/dP? stat. sys. 
[pb/GeV] 


/&±stat. 


2.5 3.3 
3.3 4.7 
4.7 15.0 


11.3 1.7 1.2 
5.47 0.67 0.60 
0.243 0.053 0.027 


0.18 ±0.03 
0.29 ±0.03 
0.27 ±0.06 


rf range 


dajdrf stat. sys. 
[pb] 


/&±stat. 


-1.3 -0.3 
-0.3 0.0 
0.0 0.3 
0.3 0.6 
0.6 1.5 


2.58 1.38 0.29 
7.98 1.55 0.88 
5.51 1.64 0.61 
7.68 1.77 0.85 
6.38 1.08 0.70 


0.15 ±0.08 
0.23 ±0.04 
0.14 ±0.04 
0.22 ±0.05 
0.21 ±0.04 


P^ u range 
[GeV] 


dajdP^ stat. sys. 
[pb/GeV] 


/ 6 ±stat. 


7 11 
11 15 

15 38 


2.82 0.33 0.31 
1.116 0.188 0.123 
0.095 0.018 0.010 


0.26 ±0.03 
0.20 ±0.03 
0.19 ±0.04 


Acft range 
[deg] 


da/dAcft stat. sys. 
[pb/deg] 


/ 6 ±stat. 


155 
155 173 
173 180 


0.0365 0.0063 0.0040 
0.462 0.056 0.051 
0.320 0.012 0.035 


0.23 ±0.04 
0.26 ±0.03 
0.11 ±0.04 



Table 7: Bin averaged differential cross sections for beauty photoproduction of dijet events 
using semi-muonic decays for x° bs < 0.75 in bins of P£, rf, P^ x , and A^j with their statistical 
and systematic uncertainties. The fit parameter f b is given including its statistical error. 
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HI Charm Dijet Muon Cross Sections 

&vis(ep —> eccX — > ejjjiX) 


P% range 
[GeV] 


da/dP? stat. sys. 
[pb/GeV] 


/ c ±stat. 


2.5 3.3 
3.3 4.7 
4.7 15.0 


49.1 4.0 5.1 
18.3 1.7 1.9 
0.854 0.126 0.088 


0.50 ±0.04 
0.47 ± 0.04 
0.31 ±0.04 


rf range 


dajdrf stat. sys. 
[pb] 


/ c ±stat. 


-1.3 -0.3 
-0.3 0.0 
0.0 0.3 
U.o U.O 
0.6 1.5 


20.3 2.5 2.1 
37.0 4.6 3.9 

42.3 4.6 4.4 

qc n /if? /in 
oo.U 4.0 4.U 

23.6 3.1 2.5 


0.44 ± 0.05 
0.45 ± 0.05 
0.48 ± 0.05 
U.4y ± U.Uo 
0.46 ± 0.05 


Pj? 1 range 
[GeV] 


da/dP^ tl stat. sys. 
[pb/GeV] 


/ c ±stat. 


7 11 
11 lo 

15 38 


11.8 0.9 1.2 
o.zz U.oy U.o4 
0.657 0.066 0.068 


0.55 ±0.04 

n on i n n A 
U.oy ± U.U4 

0.51 ±0.05 


A<^> range 
[deg] 


da/dA^ stat. sys. 
[pb/deg] 


/ c ±stat. 


155 
155 173 
173 180 


0.0820 0.0119 0.0085 
2.10 0.15 0.22 
4.06 0.35 0.43 


0.37 ±0.05 
0.52 ±0.04 
0.50 ±0.04 


x° bs range 


da/dx ^ stat. sys. 
[pb] 


/ c ±stat. 


0.0 0.4 
0.4 0.75 
0.75 1.0 


12.3 5.0 1.3 
63.5 5.8 6.6 
206.7 13.8 21.5 


0.22 ±0.09 
0.50 ±0.04 
0.51 ±0.03 



Table 8: Bin averaged differential cross sections for charm photoproduction of dijet events using 
semi-muonic decays in bins of P%, i]^, PJ? , A(ft, and x° bs with their statistical and systematic 
uncertainties. The fit parameter f c is given including its statistical error. 
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HI Charm Dijet Muon Cross Sections, a;° bs > 0.75 

o'visiep ->■ eccX ->■ ejjjxX) 


P? range 
[GeV] 


da/dP^ stat. sys. 
[pb/GeV] 


/c±stat. 


2.5 3.3 
3.3 4.7 
4.7 15.0 


31.2 3.2 3.2 
11.6 1.4 1.2 
0.368 0.119 0.038 


0.55 ±0.05 
0.49 ± 0.06 
0.20 ±0.07 


77^ range 


dajdrf stat. sys. 
[pb] 


/ c ±stat. 


-1.3 -0.3 
-0.3 0.0 
0.0 0.3 
0.3 0.6 
0.6 1.5 


14.6 2.2 1.5 
21.9 4.0 2.3 
27.0 3.4 2.8 
25.8 3.0 2.7 
9.4 1.9 1.0 


0.45 ± 0.06 
0.41 ±0.07 
0.50 ±0.06 
0.56 ±0.06 
0.39 ±0.08 


P}p range 
[GeV] 


dajdP^ stat. sys. 
[pb/GeV] 


/ c ±stat. 


7 11 
11 15 

15 38 


6.6 0.6 0.7 
3.45 0.42 0.36 
0.377 0.048 0.039 


0.57 ±0.05 
0.42 ± 0.05 
0.42 ± 0.05 


Acft range 
[deg] 


da/dAcft stat. sys. 
[pb/deg] 


/ c ±stat. 


155 
155 173 
173 180 


0.0424 0.0075 0.0045 
1.202 0.122 0.125 
3.00 0.29 0.31 


0.46 ± 0.08 
0.46 ±0.05 
0.53 ±0.05 



Table 9: Bin averaged differential cross sections for charm photoproduction of dijet events using 
semi-muonic decays for x° bs > 0.75 in bins of P£, rf, P^ a , and A<^ with their statistical and 
systematic uncertainties. The fit parameter f c is given including its statistical error. 
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HI Charm Dijet Muon Cross Sections, a;° bs < 0.75 

o'visiep ->■ eccX ->■ ejjjxX) 


P? range 
[GeV] 


da/dP^ stat. sys. 
[pb/GeV] 


/c±stat. 


2.5 3.3 
3.3 4.7 
4.7 15.0 


17.6 2.7 1.8 
6.7 1.0 0.7 
0.424 0.068 0.044 


0.43 ± 0.06 
0.45 ±0.06 
0.45 ±0.07 


77^ range 


dajdrf stat. sys. 
[pb] 


/ c ±stat. 


-1.3 -0.3 
-0.3 0.0 
0.0 0.3 
0.3 0.6 
0.6 1.5 


4.16 1.25 0.43 
12.7 2.5 1.3 
12.3 2.6 1.3 
9.9 2.7 1.0 
11.7 1.8 1.2 


0.32 ±0.10 
0.45 ± 0.08 
0.37 ±0.08 
0.32 ±0.09 
0.43 ± 0.06 


P}p range 
[GeV] 


da/dP^ tl stat. sys. 
[pb/GeV] 


/ c ±stat. 


7 11 
11 15 

15 38 


4.58 0.51 0.48 
1.216 0.375 0.126 
0.265 0.037 0.028 


0.50 ±0.05 
0.24 ±0.07 
0.65 ±0.08 


Acft range 
[deg] 


da/dAcft stat. sys. 
[pb/deg] 


/ c ±stat. 


155 
155 173 
173 180 


0.0320 0.0092 0.0047 
0.947 0.102 0.098 
0.834 0.185 0.087 


0.24 ±0.07 
0.64 ±0.07 
0.35 ±0.07 



Table 10: Bin averaged differential cross sections for charm photoproduction of dijet events 
using semi-muonic decays for i° bs < 0.75 in bins of P£, rf, P^ etl , and A<\$ with their statistical 
and systematic uncertainties. The fit parameter f c is given including its statistical error. 
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Figure 1 : The impact parameter S and Fj, el distributions for the total event sample. The data 
are compared to the PYTHIA predictions of the different quark contributions and their sum. The 
MC predictions are obtained from the fit to the total data sample. 
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Figure 2: Distributions of kinematic variables for the selected event sample. P T ea ma rac is the 
transverse momentum of the track with the highest transverse momentum in the event. The 
other notations are described in the text. The data are compared to the Pythia predictions of 
the different quark contributions and their sum. The MC predictions are obtained from the fit to 
the total data sample. 
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Figure 3: The differential cross sections for beauty photoproduction of dijet events with a muon 
as a function of P£, ff , P^ 1 , A<j^, and x° bs . The inner error bars show the statistical errors, the 
outer error bars represent the statistical and systematic errors added in quadrature. The data are 
compared with the predictions from the LO models PYTHIA, CASCADE and Herwig as well 
as to the NLO predictions of MC@NLO. The theoretical uncertainties of MC@NLO are given 
as shaded band. 27 
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Figure 4: The differential cross sections for beauty photoproduction of dijet events using semi- 
muonic decays for x° bs > 0.75 as a function of Pj,, rf, P^ etl , and Acfiti. For details see caption 
of figure |3] 
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Figure 5: The differential cross sections for beauty photoproduction of dijet events using semi- 
muonic decays for x° bs < 0.75 as a function of P^, rf, P^ etl , and A<^j. For details see caption 
of figure |3] 
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Figure 6: The differential cross sections for charm photoproduction of dijet events using semi- 
muonic decays as a function of P£, rf, -P^' 1 , A<A"> and x° bs . For details see caption of figure |3] 
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Figure 7: The differential cross sections for charm photoproduction of dijet events using semi- 
muonic decays for x° bs > 0.75 as a function of rf, -F^ 11 , and A<^. For details see caption 
of figure |3] 
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Figure 8: The differential cross sections for charm photoproduction of dijet events using semi- 
muonic decays for x° bs < 0.75 as a function of P£, rf, P^ etl , and A<^j. For details see caption 
of figure |3] 
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